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The quantum theory of conductivity of semiconductor objects, to which the quantum wells, wires and dots 
concern, is constructed. Average values of current and charge densities, induced by a weak electromagnetic 
field, are calculated. It is shown, that in both cases average current and charge densities contain two 
contributions, first of which is expressed through electric field, and second - through a spatial derivative of 
electric field. Appropriate expressions for the conductivity tensor, dependent on coordinates and applicable 
to any spatially-heterogeneous systems, are deduced. The results may be used in the theory of secondary 
light radiation from low-dimensional objects in cases of monochromatic light and light pulses. 



I. INTRODUCTION 

In connection with the increased interest to experiment 
and theory concerning light reflection and absorption by 
low-dimensional semiconductor objects - quantum wells, 
wires and dots - at light pulse excitation (see, for ex- 
ample, [ELQ) there appears a question on what kind of 
electromagnetic wave-electron interaction is more con- 
venient to use, - containing vector potential A(r, t) or 
electric field E(r, t). The work || is devoted to the same 
question with reference to calculation of a differential sec- 
tion of inelastic light scattering by infinite crystals. In 
H it is shown, that exact expressions for scattering sec- 
tions (with use sets of exact electron wave functions in a 
crystal) , obtained with use of two various kinds of electro- 
magnetic wave-electron interaction, coincide. But the set 
of exact electron wave functions in a crystal (taking into 
account precisely, for example, electron-phonon interac- 
tion) is unknown, therefore at calculation of sections the 
approached methods are used. Namely, it is taken into 
account electron-phonon interaction in the lowest order, 
not all intermediate states of an electronic system are 
considered. If to use these approximations, various kinds 
of light-electron interaction (containing A or E) will re- 
sult in distinguished results. The authors of || approve, 
that in a case of non-resonant scattering the use of inter- 
action containing E, gives the best results. 

We assume to construct the general theory of sec- 
ondary light radiation by low-dimensional semiconduc- 
tor objects, to which number quantum wells, quantum 
wires and quantum dots concern. First of all, the theory 
should describe light reflection and absorption by such 
objects, and also various kinds of light scattering (Ra- 
man scattering, Releigh scattering). The theory should 
be applicable in a case of monochromatic irradiation, as 
well as in a case of a light pulse irradiation. We are lim- 
ited to linear approximation on intensity of exciting light. 



From told above follows, that others tasks are decided, 
than in H , where light scattering in bulk crystals is con- 
sidered. Therefore the problem of an interaction choice 
(through A or E) is necessary to be solved anew. 

In the present work we shall calculate average values of 
current and charge densities induced by an electromag- 
netic field in the case of spatially-heterogeneous matter. 
To this case concern semiconductor objects of the low- 
ered dimension. Having calculated linear on electric and 
magnetic fields the contributions to average values of cur- 
rent and charge densities, we can then to determine these 
fields inside and outside of semiconductor objects, solving 
the Maxwell equations. Expressions for fields appropri- 
ate reflected and past through object light thus can be 
obtained. Such procedure, taking into account all orders 
light-electron interaction, is done in Q, where intensity 
of reflected and absorbed light at monochromatic irra- 
diation quantum well of final thickness are calculated . 
In the present work we deduce expression for average in- 
duced current density used in Q. 

The operator of interaction of the charged particles 
with electric and magnetic fields is expressed through 
vector A(r, t) and scalar <p(r, t) potentials. Therefore 
average values of the induced of current and charge den- 
sities also are expressed through these potentials. How- 
ever to use these expressions is inconvenient because of 
the contribution 

- — (0\p(r)\0)A a (r,t). 
mc 

In average of a current density , where e, m is the charge 
and mass of a particle respectively, p(r) is the operator 
of a charge density (see section IV lower). Therefore we 
shall express average values of current and charge den- 
sities through electric and magnetic fields. It is possible 
since all observable values are expressed through electric 
E(r, t) and magnetic H(r, t) fields. Our task will consist 
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in transition from expressions for physical values, con- 
taining potentials, to expressions containing fields E(r, t) 
and H(r, t). 

In the present work the case of temperature T = is 
considered. Further we assume to calculate the average 
induced current and charge densities in a case spatial - 
heterogeneous systems and final temperatures. We shall 
use some results of 

The article is organized as follows. In sections II-IV 
the statement of a task is stated, operators of current 
and charge densities and their average values on ground 
state of system are entered. In sections V -IX the task 
about expression of average values through an electric 
field and its spatial derivatives is solved. Section X is de- 
voted to exception of average values of diagonal matrix 
elements of operators r^. In section XI the general ex- 
pression for conductivity tensor of heterogeneous system 
is given. In sections XII and XIII the case of zero electric 
and constant magnetic fields is considered. 



II. STATEMENT OF A TASK. 

Let us consider s system of N particles with the mass m 
and charge e in a weak electromagnetic field characterized 
by intensities E(r, t) and H(r, t). Let us introduce the 
vector A(r, t) and scalar <p(r, t) potentials 



1 9 A ^ 
E(r,i) = — — -Vp, 
c at 

H(r,t) = V x A. 



(1) 



We shall consider fields as classical. The gage of A(r, t) 
and ip(r, t) is arbitral. For completeness of a task we shall 
consider, that the system of particles can be placed in a 
constant magnetic field H c , which may be strong (SMF). 
This field is described by the vector potential -4(r), so 



rot A(r). 



The total Hamiltonian TLtotai is written as 



n 



total 



2m * — ' \ f r 



+V(ri...r iV )+e5^(r i ,t) ) 



(2) 



Hiatal = H + U, 



(3) 



^^Ep^^-^ P* = P< - jUfo), (4) 



U = U 1 + u 2 , 

U x =-- c J d 3 rj(r)A(r,i) + J d 3 rp(r)ip(r,t), 



Un = 



2mc 



d"rp(r)A'(r,t), (5) 



where the operators of current and charge densities are 
introduced 



jM)=^ji(r,i), 

i 

g 

ji(r, t) = -^{S(r - Ti)vi + Vi5(r - n)}, v< = Pi/m, 

P( r ) = J2p^ t ^ Pi(r) =e5(r-Ti). (6) 

i 

Our task is to calculate in linear approximation on ex- 
ternal fields E(r,i) and H(r,i) the induced current and 
charge densities averaged on the system ground state. 

III. THE DEFINITION OF OPERATORS 

The charge density operator p(r) in the Schrodinger 
representation does not contain additives proportional to 
fields, but the current density operator in fields looks like 



j(r,t) + Aj(r), 



where 



A -K r ) = | E{ Av ^( r ^ r <) + <H r - ^)AvJ, (7) 



Av, = {i/K)[U,Vi] = A(r;,t), 

mc 



(8) 



[F, Q] is the commutator of operators F and Q, and hence 



Aj(r) 



-p(r)A(r,t). 



(9) 



where Vi — —ih(d/dr.i) is the generalized momentum 
operator (see, for example, ||U), ^(ri . . . r^v) is the po- 
tential energy including interaction between particles and 
the external potential. In Eq. (2) It is necessary to take 
into account non-commutativity of Vi and A{vi), A(r,-, t). 
Let us allocate in Eq. (2) energy U interaction of par- 
ticles with the electromagnetic field, having included in- 
teraction with a strong magnetic field in the main Hamil- 
tonian 



In the interaction representation we have 

p(r,t) =e* m / n p(r)e- iHt / h , 
j(r,*)=e^j( r)e -W 

Aj(r,i) = -— p(r,t)A(r,t). (10) 
mc 

Now we shall determine the current and charge density 
operators in the Heisenberg representation. In (page 
82) it is shown, that connection between the operator 
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F(t) in the interaction representation and the operator 
Fa(t) in the Heisenberg representation is expressed as 



F G (t) = S-\t)F(t)S(t), 



(11) 



where S is the S-matrix. It is determined as 

/t 
dt x U{ti) 



2 r t 



I2 



dt 2 U{t 2 ) / dt 1 U(t 1 ) + ..., (12) 



(13) 



Using Eq. (12), we find, that linear additives on poten- 
tials A(r, t) and p(r, t) to the current and charge density 
operators in the Heisenberg representation are equal 

ii Q (r,t) = Aj Q (r,i) - - j dt'U^tlU^t')], (13a) 



Pl (v,t) = --J dt'[p(r, <), (*')]• (136) 

The index 1 means the first order on potentials A(r, t) 
and p(r,t). Substituting Eqs. (9) and (5) in last expres- 
sions, we obtain 



-p(r,t)A a (r, t) 



+ - j d 3 r' j dt'\j a (v,t),j^T',t'))A (r',t') 



dV dt'\j a (r,t),p(r',t')}tp(r',t'), (14) 



U. For averaged values of current and charge densities 
we shall enter designations (0|ji(r, t)|0) and (0\pi(r,t)\0) . 
In Eqs. (14) and (15) we shall make replacement of t' by 
t" = t' — t. At averaging (0|...|0) we shall take into ac- 
count, that 

(0|F(t)|0) = (Q\e mt l h Fe- int ' h \Q) = (0|F|0>, (16) 
where F is any operator. Then we shall obtain 

(0|j la (r,t)|0) = -— (0\p(r)\0)A a (v,t) 



+ - J d 3 r' J dt'(O\\j a (r),j (r',t')]\O)Ap(r',t + t') 
-\fd\' ( dt'(0\[j a (r),p(r\t')}\OMr',t + ty, (17) 



\ Pl (r,t)\0) = -Jd 3 r'j dt'{O\[p(r,t),j (r',t')}\O) 
1 1 1 it \ 



x A (r',t + t') 

- I dV 

n 



dt>(0\[p(r),p(r',t>)]\0) 

OO 

X <p{r',t + t'). (18) 

Thus, we have obtained expressions for current and 
charge densities averaged on the ground state through 
vector and scalar potentials. But averaged values should 
be expressed through measured values - the fields E(r, t) 
and H(r, t) and their derivative. Let us make transition 
to expressions of (0|ji(r,t)|0) and (0|pi(r, i)|0) through 
fields. 



V. TIME DERIVATIVES FROM AVERAGED 
CURRENT AND CHARGE DENSITIES 



pi(r,t) = j- j dV j dt'\p{v,t),j^v'^)]A {v',t') 



— OO 
t 



£ / dV / ay[p(r,t),p(r',t')Mr',*')- (15) 



IV. AVERAGING ON THE SYSTEM GROUND 
STATE 

Let us consider a zero temperature case and average 
Eqs. (14), (15) on the ground state of our system. In all 
further calculations we shall assume, that on the indef- 
initely removed distances there are no charges and cur- 
rents, and also that on times t — > —00 E(r, t) and H(r, t) 
are equal 0, what corresponds to an adiabatic switching- 
on of these fields. In (page 84) it is shown, that at 
averaging it is necessary to use wave functions |0) of the 
ground state without taking into account the interaction 



Let us apply the following reception: let us calculate 
time derivatives from values (17) and (18): 

-(0\ji a (r,t)\0) = - {0\j la (r,t)\0) A 

+ g- t {0\j la (T,t)\0} V , 

-(0\ Pl (r,t)\0) = -(0\ Pl (r,t)\0) A 
d 



+ -(0| Pl (r,i)|0) v , 



(19) 



where indexes A and ip designate contributions of vector 
and scalar potentials, respectively, equal 

£<0|ji„(r,*)|0) x = _ J-<0|p(r)|0>^M 
at mc at 

+ ¥c J dV / dt '(°\ijc(r),3p(r',t')}\0) 

jM±t+n 

dt v ' 
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|(0| Pl (r,t)|0)A = ^ / </ V' / .//' 



x(O\[p(r),j (v',t>)]\O) 



dA p (r',t + t) 
8t 



(21) 



^(0|ji«(r,t)|0) v = -i J dV y° df 



x(0|[j Q (r,t),p(r',t')]|0}- 



0t 



(22) 



^<0|p 1 (r,t)|0> v , = -l | d 3 r'j"jt' 
dip{v',t + t') 



x(0|[p(r),p(r',O]|0)- 



(23) 



Lot us transform Eqs. (22) and (23), containing scalar 
potential <£. We use identity 

0<p(r,i + i') _ dp(r,t + t') 
dt ~~ dP ' 

then we integrate on t' in parts. We obtain 

■ J dV(0|[ jQ (r),p(r')]|0) 
x<p(r',t) 



^-(0\j la (r,t)\0) v = 



+ j: I d 3 r' I dt'(0 



Ja(r) 



,d P (v',t' 



dt' 

x<p{r',t + t'). (24) 



In the first term of the RHS Eq. (24) we execute inte- 
gration on r' and we use Eq. (6) for operators j( r ) an d 
p(v). To calculate the second term we use the continuity 
equation 



vjM ) + ^M = o, 



(25) 



which is true for operators, determined in Eq. (10) with 
taking into account the constant magnetic field H c . Fur- 
ther in this term we make integration on r' (i in parts, 
transferring derivative on the scalar potential <p(r', t+t'). 
It results in 

|(o|i-(M)|o> y = -^(o|p (r) |o)^ 



d 3 "' 



r i dt'(O\\j a (v),j (v',t')]\O) 
d<p(r',t + t') 



x- 



dr' p 



(26) 



Summing Eqs. (20) and (26) and using Eq. (1), we 
obtain 



^<0|j la (r,t)|0> = ^(0|p(r)|0>£? a (r,t) 

~\ J d 3 r' j O ^dt'{O\y a (r),j (r',t')}\O) 

xE a {r',t + t'). (27) 
We obtain completely similarly 



a(0|pi(r,t)|0) _ i 



d 3 



dt' 



dt h 
(0\[p(r),jp(r', t')\ \0)Ep{v', t + t') 



(28) 



So, we managed to express time derivatives from aver- 
aged current and charge densities through electric fields, 
having got rid of vector and scalar potentials. 



VI. AVERAGE VALUES OF CURRENT AND 
CHARGE DENSITIES EXPRESSED THROUGH 
ELECTRIC FIELD 

Integrating Eqs. (27) and (28) on time, we obtain ex- 
pressions for average current and charge densities 



( o| Mr , |o) = /_y»f 



+ C a , 



+ C. (29) 



We believe C a = C = 0, what corresponds to absence 
of induced currents and charges on times, indefinitely 
removed in the past. 

Let us enter a designation 

a(r,i) = -c / dt'E(r,t'). (30) 

J — OO 

Then with the help Eqs. (27) - (29) we obtain 

<0|ji Q (r,t)|0) - -^<0|p(r,t)|0)a Q (r,t) 
+ ^/ dV / dt'(O\\j a (r,t),j (r',t/)}\O)a (r',t'), (31) 



i 

he 



<0|pi(r,t)|0> = 
f d 3 r' f dt'(O\[p(r,t),j (r',t>)}\O)a (r',t>). (32) 

J J —oo 



Comparing expressions Eqs. (31) and (32) with Eqs. 
(17) and (18), we see, that the second differ from the first 
by absence of the scalar potential ip and replacement of 
the vector potential A(r, t) by the vector a(r, t), deter- 
mined in Eq. (30). Thus, we have managed to express 
average values of the induced current and charge densi- 
ties only through electric fields. However, without the 
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further transformations Eqs. (31) and (32) are inappli- 
cable at transition to electric field, independent on time, 
since at integration on t in Eq. (30) there appears some 
uncertainty, - the frequency u in the denominator turns 
in zero. It concerns and to the case E = 0,H^0. We 
shall return to this question in section XII. 



/ _ 



-p(r)a a (r, t). 



(38) 



It is easy to see, that Eqs. (35) and (36) may be written 



as 



j( a (T,t) = A£(r,i) - l - / dt'\j a (r,t),Ul(t% (39) 



VII. TRANSFORMATION OF EXPRESSIONS 
FOR AVERAGE CURRENT AND CHARGE 
DENSITIES 



({(T,t) = ~l dt'[p(r,t),u{(t% 



(40) 



After we managed to express average current and 
charge densities through an electric field, the task is to 
pass to expressions, obviously appropriate to electron- 
field interaction of a kind 



(33) 



which, for example, is used in M for electric field, E(t), 
independent on coordinates. 

To achieve this purpose we shall make some transfor- 
mations of Eqs. (31) and (32). We shall enter fictitious 
operators j{(r, i) and p{(r,t), for which we have 



<0|j la (r,t)|0) = (0|i/ o (r,t)|0> ) 
<0| Pl (r,i)|0) = <0|pf(r,t)|0>. 



(34) 



It follows from Eqs. (31), (32) 



-p(r,t)a a (r,t) 



+ Y c I d 3 r' I dt'[ ]a {v,t) lj0 {v\t')]ap{v',t'), (35) 



= Tcj d3r 'j_ MW,t\WJ)]af>VJ)- (36) 

Comparing the fictitious operators Eqs. (35) - (36) with 
the genuine operators Eqs. (14) - (15), we find, that for 
transition from genuine to the fictitious operators it is 
necessary to put ip(r, t) — and to replace the vector 
potential A(r, t) by a(r, t), determined in Eq. (30). 

Let us enter also a fictitious operator of interaction of 
particles with a field 



U[ = -\j d 3 rj(r)a(r,<), 



(37) 



which differs from the true operator U\, determined in 
Eq. (5) by condition <p(r, t) = and by replacement 
A(r, t) on a(r, t) .To the interaction Eq. (37) there cor- 
responds a linear on a field the additive to the current 
density operator 



what is similar to Eqs. (13a) and (13b) with replacement 
of the genuine operators by fictitious. Let us transform 
Eqs. (39) - (40) so that to remove the first term in the 
RHS Eq. (39). It is possible to rewrite the integral of a 
kind 

dt'[F(r,t),U((t')] 
from the RHSs Eqs. (39) and (40) as 

dt , [f(r,t),^ 1 (t')]-[n'. ( 41 ) 



if 



Ui = U{(t) + 



dR{t) 
dt ' 



(42) 



where F(r,t) is the operator, in case Eq. (39) equal 
j a (r,t), and in case Eq. (40) - p(r,i), and R(t) is any 
operator in the interaction representation 



R(t) = j Ht l h R Sch (r)e- iHt ' h , 



(43) 



Rsch is the operator in the Schrodinger representation. 
It is possible to show, that if 



Rsch = Rsch(?i . . . Tjv, t), 



(44) 



i.e. the operator Rsch does not contain momentums, the 
relation is carried out 



■\j a {r ) t),R(t)} = Aj Ra (r,t), 



(45) 



where 



Aj Ra (r, t) = ^e mt / h J2i^ Ra 5(r - r ( ) 



S(r - ri)Av lRa ]e' 



-iHt/n 



Avi Ra = — [Rsch, fia], 

n 



Rsch = ^[H,Rseh] H ~§f~^~' 



(46) 



(47) 



(48) 
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Also it is obvious, that under condition of Eq. (44), 

\p(T,t),R(t)] = j ut l h [ P {r),R Sch )e- im l n = 0. (49) 

Thus, it is proved, that instead of u{ (t) it is possible 
to choose any operator, determined in Eq. (42), if Rs c h 
does not contain momentums, and instead of Eqs. (39) 
and (40) to write down 

j{ a (r,t) = Aj a (r,t) - % -j dt'\j a (r,t), (50) 



p{(r,t) = -^l dt'ip^t),^')}, (51) 



where 



Aja(r) = | 5Z[A5 ia <y(r - n) + S(r - r,)Aw M ], (52) 



1 ~ 

Aw iQ = -g\Ui,r ia \. 



(53) 



Having substituted Eqs. (42) and (52) in Eqs. (50) and 
(51), we obtain 

JiaM) = -^PMKM) + l -\j a {Y,t),R{t)] 

+ ~tc J dt 'lj<*(r,t),Mr',t')}a (r',tf) 

dR(t') 



dt' 



j a (r,t), 



dt' 



(54) 



pL(m) = 

^ J d 3 r> J* dt'[p(r,t), Jf ,(r',t')}a p (r',t') 



dt' 



, \ dR(f) 



(55) 



VIII. CHOICE OF THE OPERATOR R SC h 
Let us choose the operator Rs c h as 
Rsch = - ^r i(3 a /3 (r i ,t). 

i 

Then 

fl(t) = e mt ^R Sch (r)e- mt / h 
= ~ c j d 3 rd f) (r,t)a f s(r,t) 7 

where the designation is entered 



(56) 



(57) 



d ( r ) = e^r 4 <5(r - r 4 ) = er^<5(r - r;) = rp(r). (58) 

Let us calculate the contributions in the RHSs Eqs. (54) 
and (55), containing R(t) and dR(t)/dt. It is possible to 
show, that 



1 £ 

-\j a (r,t),R(t)] = — p(r,t)a a (r,t) 
n mc 

Also it is possible to obtain the result 

dR(t) _ 1 /■ j3 _ J _,_ ^da p (r,t) 



(59) 



+ 1 /" d 3 rj {r,t)a (r,t) 



9r 7 



where 



(60) 



(61) 



Substituting Eqs. (59) and (60) in the RHS Eqs. (54) and 
(55), we find, that the terms, not containing derivative 
dap(r,t)/dt or <9a^(r, t)/dr a , are reduced. In a result we 
have 



ria = — <w>*)- 



TlC 



3„' / A+l 



dt'\j a (r,t),d {r',t')] 



■ , .. dapjr, t) 
m<- ' dr a 

, , da (r',t') 



dt' 



-i- c J rfV j dt'\j a (r, t), W, t')] f) , (62) 



he 



dt , [p(r,f),d /J (r , ,t')] 



pf(r,t) = 



Of 



~ J d 3 r' dt'[p{v, t), Y 01 (r>, t')] ^/ *' } . (63) 

Averaging the fictitious operators Eqs. (62) and (63) 
on the ground state, we shall obtain required by us ex- 
pressions for averaged values of the induced current and 
charge densities. 

Taking into account the definition Eq. (30), we find, 
that terms in the RHS Eqs. (62) and (63) are divided 
on two categories: to the first concern what contain an 
electric field, to the second - what contain derivative from 
this field on coordinates. Therefore it is convenient to 
write down average values from Eqs. (62) and (63) as 

(0|ji Q (r,t)|0) = (0Ui Q (r,t)|0>B + (0Uia(r,t)|0)s B/8l ./ ) 

(64) 
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(0|pi(r,t)|0) = (0|pi(r,t)|0) B + (0\pi(r,i)\0} dE/dr , (65) 
where 

(0\ Jla (r,t)\0) E = 

l -(d\' f dt'muv^.dpiv'xmEpiv^t 1 ), (66) 



{0\ Pl (r,t)\0) E = -j dV 
x / dt'(O\[p(r : t),d (r',t')}\Q)E p (r',t'), (67) 



(0\ 3la (r,t)\0) dE/dr = --(0\d p (r)\0) [ 



+ - I d 3 r' I dt'(O\\j a {r,t),Y 01 (r',t' 



x / dt 



„dEp{r>X) 



(68) 



(0\pi(r,t)\0) eE/dr = ~Jd 3 r'J dt' 

/*' f)E (r' t"\ 
dt" - ^/ ' . (69) 



So, the task of expression of averaged on the ground state 
current and charge densities in linear approximation on 
electric field and its derivatives on coordinates is solved. 



IX. THE TRANSFORMED FICTITIOUS 
OPERATOR OF INTERACTION 

For completeness of a picture we shall determine a kind 
of the fictitious operator of interaction U\ = u{ + Rsch- 
Using Eqs. (37) and (60) we obtain in the Schrodinger 
representation 



U\ — Uie + U ldE / dr , 



(70) 



dE^t') 



dr. 



*7 



Or 



27 I ' 2/3 



dEpjn,?) dEpjr^t') 

+ n/3 — ^ 1 ^ v 



dr v 



[ dt'^ivi 

J - OO i I 



<9r 47 " ll 
dEpjr^t') 



dr 



Ti/3 



dEpjr^t') 
r i0 ^- v n }. (,4) 



Thus, U\ E contains an electric field, and Uig E /g r con- 
tains derivatives from electric field on coordinates and 
an integral on time. In the case, when terms, containing 
derivatives from electric field on coordinates, on what of 
reasons give the small contribution □, it is possible to use 
Eq. (73). It coincides with the formula for interaction of 
charged particles in |p|,p|,p0| . Let us notice, that Eq. (70) 
for fictitious interaction U\ may be written compactly 



< ^ I da p {r t ,t) dap(ri,t) 
2c ^ I dt np + 7-1/3 dt 



(75) 



where dap(ri,t)/dt is the full time derivative: 



dap{ri,t) da {ri,t) i, , . . 

— = + £ [«><*(*«>*)]. (76) 



and 



dap(ri,t) 



-[W, fl ,(r„t)] = -|, 7 — 
dap{Yi,t) ^ 



(77) 



The following linear on a field additive to speed corre- 
sponds to interaction Eq. (70) 



e dap(r it t) 
Av ia = — r i/3 — . 

mc OTia 



(78) 



what corresponds, according to Eq. (52), to an additive 
to the current density 



<^rdp(r)Ep{r,t), 



UidE/dr = - \ d J rYa 7 (r) / dt 



dr~< 



(71) 



(72) 



Having executed integration on r in the RHS Eqs. (71) 
and (72), we obtain 



Uie = -e} j r i pEp{Y h t), 



(73) 



Aj Q (r) = ^(r)^M. 
mc ar a 



(79) 



Passing to the interaction representation, we obtain the 
first term in the RHS Eq. (62). 



'Strictly speaking, it follows from the Maxwell equations, 
that if the electric field depends on time, it also depends on 
coordinate, i. e. derivatives of field components on coordi- 
nates are distinct from zero. 
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X. EXCEPTION OF DIAGONAL ELEMENTS OF 
OPERATORS OF PARTICLE COORDINATES 

Let us return to Eqs. (64) - (69) for average induced 
current and charge densities at T = 0, obtained in sec- 
tion VIII. Taking into account definition (58) of operators 
d(r), and also definition (61) of operators Yp 7 (r), which 
may be rewritten as 



(80) 



_ mt/n dap(rj,t) ^ mt/h 



dt 



(85) 



and the derivative dap (r, , t) j dt is determined in Eq. (76) . 
Having executed in the first term of the RHS Eq. (83) 
integration on t' and having calculated the commutator 
\j a (r),ap(ri,t)], we obtain the result 



x ia p{r,t) = 0. 



(86) 



we see, that these operators contain coordinates of par- 
ticles. But averaged values (0|ji(r, t)\0) and (0|pi(r, t)\0) 
do not owe to depend on an initial point of readout of 
coordinates r,. Let's transform expressions (64) - (69) so 
that last property became obvious. Let us represent a 
vector rj as two parts 



r t =v t + (0|r l |0). 



(81) 



It is obvious, that matrix elements of the operator fj, 
as diagonal, and not diagonal, do not vary at change 
points of readout of coordinates iv We shall show, 
that in expressions for average values (0|ji Q (r, t)|0) and 
(0|pi(r, £)|0) the operators may be replaced by f j. Let's 
write down average (0\jia(r, t)\0) as 

(0|ji a (r, t)\0) = (0\j la (r,t)\0) + ^(0|r l/3 |0)a; M/3 (r, t), 

i 

(82) 



where (0|ji Q (r, i)|0) - is the contribution of the operators 

x ia p(r,t) = l -J d 3 r' J dt'(0\{j a (r,t), Pi (v',t')}\0) 
x^,O + ^(0|«W|0>^ 



- 4- d 



he 

da/sir', t>) 
dr' 



dt'{0\[j a (r,t), Jl7 (r',t')}\0) 



(83) 



Having executed in the first and third terms of the RHS 
Eq. (82) integration on r', we obtain 



x ( 



ie f* 

x ia p(r,t) = -— dt' 
dap(r t (t'),t') 



e 



dt 1 



(84) 



(0|ji a (r,t)|0) = (0|j lQ (r,i)|0). 
It is similarly possible to obtain 



(87) 
(88) 



(0|pi(r,t)|0) = <0|pi(r,t)|0). 
Let us represent Eqs. (87) and (88) as two parts 
(0|ji a (r,t)|0) = (0|ji a (r,t)|0)/ + (0|ji a (r, t)\0)n, (89) 

<0|pi(r,t)|0> = (0|pi(r,i)|0)j + (0|pi(r,i)|0)//, (90) 
where 

(0|ii a (r,t)|0>j = 
\j^rij dt'(O\[j a (r,t),d0(r',t')]\O)E f3 (r , ,t'), (91) 



(0\j la (r,t)\0) n = £< |^(r)|0)^^ 

lit C CJTq 



ch 



f d 3 r' f dt'{Q\\j a {v,t)Y M {v',t')\\Q) 



dr' 



(92) 



(0|pi(r,t)|0) 7 = 
Ijdtr'f df(0|Kr,t),^(r',f)]|0>^(r',O. (93) 



{0\ Pl { r ,t)\0) II = ~ J d?r' f dt' 

dap{v',t') 



and 



where the designation for the full time derivative is used 



x(0|[p(r,t),^ 7 (r',i')]|0)- 



d(r) = e^r iP (r), 



(94) 



(95) 



(96) 
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Results (89) - (94) are basic in the present work. Let's 
emphasize, that splittings Eq. (89) and Eq. (90) average 
sizes on two parts do not coincide with splittings Eq. (64) 
and Eq. (65). 

The contributions with an index I we shall name basic, 
as they do not disappear in case of an electrical field in- 
dependent from coordinates r. Contributions // contain 
derivative from an electrical field on coordinates. 



XI. THE CONDUCTIVITY TENSOR, 
DEPENDENT ON COORDINATES 

Let us consider at first only basic part of the induced 
current density, designated in Eq. (90) by index /. We 
shall write down Eq. (91) as 

(0|ji a (r,t)|0)j = 
= f d 3 r f dt'l af3 (v,t,r',t')Ep(r',t'), (97) 

J J — OG 

where 

l a p(r, t, r', t') = % -&{t - i')(0|[ja(r, t),dp(r', t')}\0). (98) 

Let us enter a tensor 

a Ia p{r', t'\r, t) = l af3 {r, t,r-r',t- t'). (99) 

Designation Eq. (99) with a partition we have borrowed 
from . Then it is possible to rewrite Eq. (97) as 

(0|ji a (r,t)|0)j = 
= fd 3 r'[ dt'a Ia p(T',t'\r,t)E p {r-r',t-t'), (100) 

J J OG 

where 



a Ia/3 (v',t'\v,t) = -e(t')(O\\j a (r,t),d (r-r',t-t')]\O). 

(101) 

It is visible from Eq. (101), that the tensor cr/ Q , i g(r', t'\r, t) 
does not depend on t\ 

Now we shall make the Fourier transformation. Let us 
write down electric field as 



£ Q (r,i) = £W(r,i) + £r J (r,i), 



(-), 



(102) 



where 



E^(r,t) = J d 3 k J™ dcoE a (k,u;)E 



Ei-Hv,t) = (Ei+\r,t)Y, 



(103) 



(104) 



E a (k,u)= ( d 3 r [ dtE a {r,t)E 

J J — OO 



— ikx-\-iujt 



(105) 



Let us enter a Fourier-image of the tensor cr/ Q ^(r', i'|r, 0) 
on variables r' , t' 



<7 Za/ 3(k,u/|r) = f d 3 r' ( ai a p(r',t'\r,0) 

J J — OO 



(106) 



Then 



(0|jia(r,t)|0) 7 = <0|ji a (r,t)|0)y 



(+) 



(107) 



(0b' la (r,t)|0)( +) = 
^4 Jd 3 k J d W a /Q/3 (k,^|r)£; / 3(k, W )e lkr - i " t , (108) 



(0|j la (r,t)|0)^ ) = ((0|j la (r,t)|0)f0*- 
Substituting Eq. (101) in Eq. (106), we obtain 



(+h 



(109) 



f\ —i\cr' -\-icjt' 



a Ia0 {k,Lj\r) = - I d 6 r' I Q(t')e 

x(0\\j a (r),dp(r-r',-t')]\0}. (110) 

By similar way we find the contribution in conductivity 
with an index //. Finally we obtain 

(Q|j a (r,i)|0>« - 

w J d3k So '^^^z^^y^^' ( in ) 



(2 

cr Q /3(k, w\r) = ai a p(k, ui\r) + a Ha/3 (k, cu\r). (112) 
The value erj Q/ 3(k, ui\r) is defined in Eq. (110), 

a 7/Q/3 (k, W |r) = ^(0|J /3 (r)|0) 
moj 

[ d \' [°° dt'Q(t')e- iW+iut ' 
muj J J_ 00 

x<0|[7a(r),y /j7 (r-r',-t , )]|0). (113) 

The conductivity tensor <7 a p(\i, cu\r) does not depend 
from coordinates r only in case of spatially - homogeneous 
system. At consideration of semiconductor objects of the 
lowered dimension dependence of the conductivity tensor 
from r is rather essential. Dependence of the conductivity 
tensor of the low dimensional semiconductor objects on 
r is rather essential. 

In all our previous works (see, for example, 
12|jl3 14]) Eq. (91)was used for calculations of an in- 



duced density of a current. 
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XII. TRANSITION TO EXPRESSIONS 
CONTAINING MAGNETIC FIELD 

Till now we left behind frameworks of our consider- 
ation a case, when electric field E does not depend on 
time, in particular the case E = 0, H(r,i) = const. 
To consider the last case we shall transform expressions 
for the average induced current and charge densities ob- 
tained at the end of section VIII, having entered in them 
magnetic field H(r, t). For this purpose each of values 
< 0|ji Q (r,i)|0 > aE /dr and < 0|pi(r,t)|0 > aE /dr, deter- 
mined, respectively, in Ens. (68) and (69), we shall break 
on two parts as follows 111 



(-) 



(0\ji a (r,t)\0) dE/dr = 
(0|iia(r,t)|0)W + (0|j la (r,t)|0)(-), 

(0\ P i(r,t)\0) dE/9r = 
,(0| Pl (r,i)|0)( + ) + (0| Pl (r,t)|0)(-), 



(114) 



(115) 



where 



<0|j la (r,i)|0> 



(±) _ 



;(0|^(r)|0) 



2mc 

/ dap{r,t) da a (r,t) 
V dr a drp 

dt'(Q\[j a (r,t),Y^(r',t')}\0) 
( da {v\t') ± d ai {v',t')\ 



dr' 



(116) 



(0| Pl (r,i)|0) ± = -^- / ,h' I ,//' 



— oo 

/ 4-1 



x (0\[p(r,t),Y^(r',t 
± 



cMr',t') ± da 7 (r',t') 
dr' dr'„ 



(117) 



At first we shall consider contributions with index (-). 
Let us return to vector A(r,t) and scalar <p(r,t) poten- 
tials. Taking into account definition Eq. (30) for the 
vector a(r, t) and first of the formulas Eq. (1), we obtain 

a(r,i) = A(r,t) + c f dt'd(f(r',t')/dr. (118) 

J — oo 

Having substituted Eq. (119) in the RHSs of expressions 
(0|ji a (r,t)|0)(-) and (0|pi(r,t)|0)(->, we find, that con- 
tributions from scalar potential ip become zero, and 



^Here top indexes (+) and (-) have no any relation to the 
same indexes in section XL 



0bia(r,t)|0) 

2mc 

(dA p {r,t) dA a (r,t) 



(0|dfl(r)|0> 



i 



Or 



, )//; , dr' / dt'{0\\j a (r,t),Y^(r',t')]\0) 
'dAp{r',t') dA y (r',t') 



dr' 



dr' 



(119) 



x (O\[p(r,t),Y 01 (r\t')}\O) 
fdAp(r',t') dA. ( {r',t') 



dr' y 



dr' (j 



(120) 



Taking into account, that H(r, t) = rotA(r,t) (second 
equality from Eq. (1)), we obtain easily 



<0|ji a (r,i)|0) 



(-) _ 



2mc 



(H(r,t) xr) Q (0|p(r)|0) 



+ Wc J^J_jm(r',t'),r % 

x(0|[j Q (r,t), J/3 (r',i')]|0), (121) 

(0\pi(r,t)\0)- = ^-Jdr'£jt> 
x(H(r',i') x r>) p (0\{p(r,t), 3 p(r',t>)}\0). (122) 



Thus, the values (0\j la (r,t)\0)^ and (0|pi(r, t)\0)^ are 
expressed through magnetic field. 

Now we shall transform Eqs. (117) and (118) for the 
contributions with an index (+). We work under the 
following circuit. Let us consider, for example, Eq. (117). 
In the RHS part 

W r = ^i-y(r') (123) 

enters; for j'-y(r) it is easy to deduce a ratio 

j 7 (r) = dy(r) + dY l5 {r)/dr s . (124) 

Let us substitute Eqs. (124) and (125) in the second 
term of the RHS Eq. (117), which in result breaks up 
on two contributions - occurring from dY 7 g(r')/dr' s and 
from d~f (r) . In the first contribution we integrate in parts 
on r' s , in second - on variable t' also in parts. It results 
in 

<0|j la (r,t)|0>(+) = -(0|j la (r,t)|0>W 
-— <0|p(r)|0)r p r 7 



+ — I dr'r'^ I dt'(0\\j a (r,t),j s (r',t' 
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times 



dr' 7 dr' 5 he 



dr'r' a / dt' 



'0 

12„ + l\ 



x (0|ba(r, t), d, (r', t'W) 9 a d fgf ] - (125) 



7 



Using a ratio 

-c- 1 dap(r,t)/dt = E p (r,t), 
we obtain finally from Eq. (125): 



2mc 
/ 



dr a dr 7 



+ 2he J dr ' r '^ y_ oo ^'(0|L7a(r,i),^(r',t')]|0> 

^(r^) i /■ " 
9^ 2ft 7 d r ^ 

x (0\\j a (r, t),p(r',t')]\0) SLJ L , (126) 



dr> 7 



and, similarly, 

<0|p(r,<)|0)(+) = ^-|drV^; ^' df 

x(0|[p(r,t), j5 (r',t')]|0)^^ 

X (0\[p(v,t),p(v',t')]\0)^^l. (127) 

Let us notice, that unlike Eqs. (68) and (69), Eqs. (127) 
and (128) contain only second derivative from a vector 
a(r,t). The obtained results can be written down in the 
symmetric form: 



(0|j la (r,f)|0>(+) 



4mc 



{0\p(T)\Q)r p r n 



d i 
x g^r {dap (r , t) /<9r 7 + da y (r, t) /drp ) + — 



x dr'r'prU dt'(0\\j a (r,t),n^(r' ,t)}\0), (128) 



(0|p(r,t)|0)W = ^ 
xfdr'r'^f dt'{Q\[p(r,t),Slp y (r',t')]\0) (129) 



XIII. CONSTANT MAGNETIC FIELD 

Let us consider a case of a constant in time and space 
magnetic fields H = const, E = 0. Let us remind, that 
we have included constant magnetic field in unperturbed 
Hamiltonian H (see section II). But in the present section 
we do not include a field in our Hamiltonian and consider 
it so weak, that it is possible to be limited by linear on 
field contributions to induced current and charge densi- 
ties. As E = 0, contributions to the average induced by E 
densities are equal ( See Eqs. (66) and (67)). The con- 
tributions with indexes dE/dr are broken on two parts, 
which we have designated by the top indexes (+) and 
(-). Parts with indexes (+) are equal 0. It is easy to be 
convinced of it if to choose vector and scalar potentials, 
for example, in gage || 



A(r) = (1/2) [r x H], ip = 



(131) 



and to use Eqs. (118), (126) and (127). 

So, there are only contributions with indexes (-), deter- 
mined in Eqs. (122) and (123). Having put H(r,t) = H 
and having made replacement variable t' by t" = t' — t, 
we obtain time-independent results: 

(0\ji a (r,t)\0) H = -^(H x r) Q (0|p(r)|0) + JL 
x [ dv 1 f dt"(Hxr')p(0\\j a (r),jp(r',t"W), (132) 



<0|MM)|0) ff = — 
x I dv' I dt"{KxT')p(0\[p(v) l]p {v\t"m. (133) 



The indexes H in the LHSs mean H = const. E = 0. 

In obtained expressions we shall pass from operators 
Ti to the operators Fj, determined in Eq. (81). By a way, 
similar to a way stated in section X, we obtain 

<0|ji«(r,t)|0>H = -^-(H x <0|d(r)|0» Q 
Zmc 



2kc 



dt ,J £{(0\\j a (r,-t'),(U x f i)pv ]\O), (134) 



•\pi(T,t)\0) H 



is 
2hc 



where the symmetric tensor 
d 

tyg 7 (r,i) = j s (r,t) — (dap(r,t)/dr y + da^(r,t)/drp) 

+ p(r, t)^- t (dap(r, *)/0r 7 + da 7 {r, t)/dr ), (130) 



x / ^'^{(Olb^-tO^Hxr^^jlO). (135) 

J-oo - 

Let us note that the relations are performed 
div(0\h{r,t)\0) H = 0, / (0|ji(r,t)|0) H dr = 0. 
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XIV. CONCLUSION 

Let us consider the basic results. Since operator of in- 
teraction of charged particles with electromagnetic field 
is expressed through potentials A(r, t) and ip(r, t), but 
not through electric E(r,t) and magnetic H(r, t) fields, 
initial expressions for average values of induced current 
and charge densities also are expressed through poten- 
tials (see Eqs. (17) and (18)). But Eq. (17) for a current 
density is inconvenient, since it contains the contribution 
proportional to average (0|p(r, t) |0) of a charge density, 
and this contribution is not completely small. 

Therefore we have put a task to express induced den- 
sities through fields E(r, t) and H(r,<), that is certainly 
feasible, since the average current and charge densities 
are observable values. 

In result Eqs. (31) and (32) were obtained, containing 
only electric field. However, Eq. (31) has the same lack, 
as Eq. (17), - it contains the contribution proportional 
(0|p(r,t)|0). 

If to use an approximation, in which the electric field 
E(r, t) does not depend on coordinates, but only on time, 
the contribution proportional (0|p(r, t) |0) , may be re- 
moved, if to enter in expression for average current den- 
sity the operators of «th particle coordinates , having 
taken a ratio = (i/K)[H, r-j]. This reception in essence 
is used in || , where the theory of light scattering in bulk 
crystals is discussed. In the appendix A it is shown, how 
the same reception allows to exclude the charged parti- 
cles concentration from the formula for conductivity of 
intrinsic bulk semiconductors at T = 0. 

The same result turns out, if in approximation 
E(r, t) ~ E(t) to write down interaction of particles with 
a field as — e J2i r i(3E^(t), as it is made in 

However, the task becomes complicated if to take into 
account dependence of an electric field on coordinates. 
We have put to ourselves by the purpose to obtain such 
expression for induced current density, which would pass 
in the Kubo-formula at transition from E(r,t) to E(i), 
and similar expressions for induced charge density. For a 
case T = this task is solved in sections VII - IX. 

Essence of used reception consists in the following. The 
operator of interaction of particles with a field, expressed 
through electric and magnetic fields, does not exist. But 
we enter a fictitious operator of interaction, which results 
into correct results for average on the ground state of sys- 
tem values of the induced current and charge densities. 
This fictitious operator ll{ Eq. (37) is expressed only 
through an electric field E(r, t), in an integral on time. 
Then such transformation of interactions Ll{ — > U\ is 
made, which does not change values of average densities, 
but eliminate the contribution containing (0\p(r,t)\0) 
from expression for average induced current density. In 
result Eqs (64) - (69) for (0|ji(r, t)|0) and <0|pi(r, t)|0) 
are obtained, in which the basic contributions contain an 



electric field, and additional - the derivatives from E(r, t) 
on coordinates. The transformed fictitious interaction U\ 
also is divided on basic and additional parts. The first is 
equal — e ^ r^iJg^, i), the second contains derivative 
dE a (r,t)/dr/3 in an integral on time. 

In section X Eqs. (64) - (69) for average densities of 
induced currents and charges are transformed so, that 
obviously do not depend on an origin of coordinates of 
particles. It is the basic result of the present work. In sec- 
tion X from Eqs. (64) - (69) for average induced currents 
and charges densities the diagonal elements (0 1 1 0) are 
excluded. The results Eqs. (90) - (96), containing only 
non-diagonal elements of the operators r^, are obtained 
. This is the main result of the present work. In section 
XI the expressions for the conductivity tensor cr(k, u\r) 
are obtained, dependent from coordinates, in a case of 
spatially-heterogeneous systems. 

Further in section XII from expressions for (0|ji (r, t) 0} 
and (0|pi(r, t)\0), containing only electric field and 
derivative from i) on coordinates, we pass to expressions 
containing also magnetic field H(r, t). It is necessary for 
reception of results in case H = const, E = 0, which is 
considered in section XIII. 

In the Appendix B the expression for the operator of 
total acceleration of system of the charged particles is ob- 
tained. The acceleration is caused by an external weak 
electromagnetic field. It is shown, that averaged on to 
the basic state of system the total acceleration can be 
expressed through electrical and magnetic fields. In the 
case of free particles the obtained result passes to a cor- 
rect limit containing a force caused by an electric field, 
and the Lorentz force. Further it is supposed to use ob- 
tained results for construction of the general theory of 
secondary radiation from low-dimensional semiconductor 
objects. 

The work has obtained the partial financial sup- 
port from the Russian Fund of basic researches (00- 
02-16904), Program MNTK " Physics of semiconductor 
nano-structures" and the Federal Program " Integration" . 

APPENDIX A 
Three expressions for the conductivity tensor in 
a case of a spatially - homogeneous matter and 
an electric field independent on coordinates. 

Proceeding from Eq. (31) for the average induced cur- 
rent density and taking into account a following connec- 
tion between the Fourier-components 

a(k,w) = (ic/o;)E(k,a;), (Al) 

following from definition Eq. (30), we obtain 

ie Iff 00 
<7^(k,w|r) = (0|p(r)|0)* a/s + — / dr' / dt' 

Xe -*r'+^ e(t / )(0 | iQ(r)!j>(r _ r , _^ m _ {A2) 
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In a case of spatially - homogeneous system the tensor 
<r aj g(k, u\r) does not depend on r. In the first term we 
use a ratio 



(0|p(r)|0) = en, 



(A3) 



where R a — J2i r ia- Having substituted (A10) in (A8), 
we obtain 

ie 2 n e 2 f°° 

o a p(u) = 5 a p + — — / dtQ(t) 

mu huV J_ 00 



where n is the concentration of charged particles, in the 
second term we integrate on r and r', also we shall divide 
the result on the normalized volume Vq. We obtain 



xe 



iu)t 



dt 



(0\[R a (t),Vp}\0). 



(All) 



ie 2 n , 



a a p(k, u) = 5 a p 



mu 4m 2 huVo 



/ dtO(t)e 



iut 



Let us execute integration on t in parts. At t —> oo 
function Q(t)e iut -^(0\[R a (t),Vp}\0) -> at replacement 
u by u + iS, 5 — > 0. Let us take into account also, that 
dG(t)/dt = S(t) and 



> /{] ^ e iHt/h( e -ter ipia 



(l/V )[R a ,Vp] = (ihn/m)S a0 , 



(A12) 



then we shall obtain 

oo 



where p 4 is the momentum operator determined in Eq. 
(4) with the account of a constant strong magnetic field. 

In case of electric field E(i), independent on coordi- 
nates, we enter the frequency representation 

/oo 
dtE a (t)exp(iojt), (A5) 
-oo 



and 



E a (k,co) = (2ir) 3 5(k)E a (uj). 



(A6). 



Let us enter also conductivity a a p(u\r), included in def- 
inition of an average current density 



(0| Jla (r,t)|0)+ = 



1 

2^ 



dwa a j3 (to | r ) E a (to ) e 



(A7) 



where the index h means spatially-homogeneous electric 
field. 

Using section XI, it is easy to show, that 



a a/3 (uj\r) = cr Q/3 (k = 0,w|r). 



(AS) 



From (A. 4) and (A. 8) in a case of spatially - homogeneous 
matter and the field E(t) we obtain 



ie 2 n 



muj 



5 a f 



e 2 I 00 

— J Jte(t)e^(0\[V a (t),V,}\), (A9) 



+ 



where V a = (l/m)J2iPia is the operator of total speed 
of charged particles. (A9) is the first formula for the 
conductivity tensor. For reception two others we use a 
ratio 



V a = (i/h)[H,R a ], 



(AW) 



where the designation is entered 
D = eR. 



(AU). 



Eq. (A13) is the second expression for the conductivity 
tensor. Let us emphasize, that at transition from (A9) 
to (A13) there was a reduction of the first term from the 
RHS (A9), containing concentration n. 

At transition to the third expression (A10) for the 
operator Vp from the RHS (A13) is used. Further we 
work under the same circuit, as at transition from (A9) 
to (A13). Taking into account, that the commutator 
[D a ,Dp] = 0, we obtain the third formula 



C 7 M 



dte(t)e^(0|[^ Q (t),^]|0). (.415). 



It is convenient to apply (A9), (A13) and (A15) to 
different systems. 

For free particles, when V(y\...vm) =0, H c = 0, 
the velocity operator v commutes with the Hamiltonian 
'Hfree = m J2i v i • Let us use (A 9), taking into account, 
that v a (t) — v a , [v a ,vp] — 0, and we obtain the known 
result 



ie 2 n 
mu 



(A16) 



In a case, when exited states are separated from the 
ground state by an energy gap , i. e. the energy of these 
states E n = fiu n ^ 0, it is more convenient to use (A 13) 
or (A 15). Using exact wave functions |n > of the exited 
states and having calculated integral on time, we shall 
express (A 15) through matrix elements of the operator 



D: 



n. 



< 



0\D a \n >< n\Dp\0 > 

u - u n 



13 



< 0\Df,\n >< n\D a \0 > 



LO + LO ri 



(A.17) 



whence follows, that at to — > a a p(ui) — > for systems 
with an energy gap (superconductivity is not considered) . 
To number of systems with an energy gap the bulk semi- 
conductors without impurity and defects concern. 

The conductivity tensor <j a p{uo) is connected with the 
dielectric susceptibility tensor e a p(ui) by ratio 

e a p(u) = 5 a p + 4irx a f3{u), 

so from (A.17) we obtain 

. 1 / <0\D a \n><n\D p \0> 



+ 



< 0\Dp\n >< n\D a \Q > 



(A.18) 



LU n +LU 

Using the ratio between matrix elements 

< 0\D a \n >= (ie/uj n ) < 0\v a \n >, 

< n\D a \0 >= -(ie/oj n ) < n\v a \0 >, 
from (And 18) we find 

< 0\v a \n >< n\v \O > 



LO n - LO 



+ 



< 0\vp\n X n\v a \0 > 



■ LO 



(A19) 



If w < w„, Xa/3 is real and also does not depend on 
frequency lo, but at lo ~ u) n this dependence becomes 
strong and there appears an imaginary part, distinct from 
zero, of the tensor Xup, determining the resonant light 
absorption on frequencies to ~ LO n . For calculation of 
an imaginary part it is necessary to replace frequency 
lo by lo + ry/2, 7 — > 0, when the field is switched on 
adiabatically, or to replace frequency w n by ui n — i~//2, i. 
c. to take into account a final life time of a system in a 
state n. 

APPENDIX B 
Acceleration of particles system 

To obtain of the operator Ji(t) of the induced current 
of system of particles we use Eq. (14), containing vector 
and scalar potentials. Having executed integration on r 
and r', we obtain 



J la (t) = --V4„(ri(t),t)| 
mc £ -~' 



-n 1 dt 



(SI) 



where the designation is used 

A(Ti(t),t) = e m */ fi A(r,,i)e- 4 ' Ht /' i . 

Differentiating (Bl) on time and having divided on e, 
we shall obtain the operator Wi(t) of the induced total 
acceleration: 



W la (t) = --^Vifr.fi),!) 
mc at /L — 1 



tfi(i),5^« ia (t) 



+ 1 f dt'iU^W^t)}, (52) 

where W a (t) = Y,i w ia(t), w ia (t) = dv ia (t)/dt. The 
decoding of first two terms gives 

W m- e \^ dAa(ri(t),t) e ^ dip{vj(t), t) 

la[) mc^ dt m^ dr ia (t) 



+ £ {«*(*)( 



dA {Yi{t),t) dA a (ri(t),t) 



dr ia (t) 



dr i/3 (t) 



f dA p ( ri (t),t) _ dA a ( ri (t),t) \ 
+ \ dr ia (t) dr ifi {t) ) Via[t) 



+ 



(i/ft) f dt'lU^t')^^)}, (53) 

J —00 



where the designations are used: 

dAa{Ti(t),t) _ iHt/h dA a (ri,t) mm 
8t dt 6 ' 



d<p(rj(t),t) 
dv ia (t) 



iHt/h ( r i > *) r -iHt/H 

dr in 



H c is the constant strong magnetic field included in the 
basic Hamiltonian H. At transition from (B2) to (B3) 
the ratio is used v, x v, = (ihe/m 2 c)H C7 and the LHS 
7^ because of non-commutativity of various projections 
of the velocity, for example 



[vi x Wi] z = [v ix ,v iy ] = (ihe/m 2 c)H cz . 



(B4) 



The operator w, (t) of acceleration ith particle is equal 



w ia = U iHt ' h [H,v ia ]e-™/' i = 



14 



1 cl m/h dv (ru---rN) c ^ m/h 

dr ia 



m 



+ 



e(vi(t) x H c ) c 



toc 



(55) 



With the help Eq. (1) it is easy to see, that the expression 
in braces in (B3) is equal 



mt/h 



{(v, x H(r t ,t)) a - (H(r 4 ,t) x Vj) a }e 
Therefore (B3) will be transformed to 

Wi(t) = -^£ E M*)>*) 



-iHt/n 



+E l3 {r i (t'),t + t')v jl3 (t')} 



(58) 



This expression contains only electric fields. On the other 
hand, having averaged the operator (B6), we obtain 



(0\W la (t)\0) = -y2(0\E a (r u t)\0) 



+ ^;5Z(0|(v 4 x U(n,t)) a - (H(ri,t) x v;)a|0) 



-(e/m C ) 2 ^(0|(A(r,,t)xH c ) Q |0)+C Q (t), 



(59) 



+^- c E< v *w x H ( r « (*)>*) - h *m^ *) x v * (*)> 



(£) 2 EA(rW,t)xH c 



where 



c a (t) = (i/h) f dt'<o|[t/i(t'),w' Q (t)]|o). (BIO) 

J — CO 

Let us transform (BIO). Integrating in parts, we obtain 
C a (t) = Ci(t) + C 2 a (t), (511) 



dtTOtO.Wti)]. 



(B6) 



It is obvious, that the second term corresponds to the 
Lorentz force, at which the non-commutativity of the 
operators and H(r i7 1) is taken into account. The third 
term caused by H c , contains the additive to speed Avj, 
determined in Eq. (8) and induced by a weak electro- 
magnetic field. 

In the case of free particles 



H c = 0, V{t u ...t n ) = 0, W = 0, 



(B7) 



and in the RHS (B6) only two first terms containing weak 
electric and magnetic fields are kept . 

But if particles are not free, the operator (B6) fails 
to be expressed only through fields, since the last two 
terms contain vector and scalar potentials. Average value 
(0|Wi(t)|0) of the induced acceleration should be ex- 
pressed only through fields, that we now shall prove. For 
this purpose we shall calculate the value (0|Wi(i)|0) in 
another way, and then we shall check up, that both ways 
give coinciding results. Let us use Eq. (27) and integrate 
its both parts on r. We obtain 



<0|W lQ (t)|0) = -V(0|5 Q (r 2 ,i)|0) 

TO 



C 1 a (t) = -(i/h)(0 



cl(t) 



i d f 
Udl 



f dt' 

J -co 



(512) 



) . (513) 



In the RHS (B12) we shall substitute Eq. (5) for Ui and 
we shall calculate the commutator. We obtain 



d<p(r u t) 



dn, 



Kmc 



J 5>|(A( ri ,t) xH c ) c 



2toc^ \ 



dA p {Ti,t) , dAp(n,t) ^ 



dn, 



0). (514) 



In (B13) in integral we pass to t" =t — t' and we break 
the integral on two parts, first of which contains vector, 
and second - scalar potentials: 



cl{t) = c\ a {t) + cl a (t), 



(515) 



~2h. 



dt' ( 



c 2 Aa (t) 




E" 
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dA p {v t {t'),t + t') 



dt 



dAp( ri (t'),t + f) ^ , 
H ^ «»/H* j 



dt 



df 



x ( 



9^(0, t + o 



(516) 



(517) 



We shall leave (B16) without changes, and in (B17) we 
use a ratio 

dip{ri(t'),t + t')/dt = 
= exp(mt' /h)(d<p(Ti,t + t')/dt) exp(-iHt'/h) = 



Summing (B14), (B16) and (B20), we obtain finally 



C a (t) 



2mc 



Vi0 



dAp{r u t) dA fi {v t ,t) 



dr it 



+ 



On, 



+ (e/m C ) 2 ^(0|(A(r,,t)xH c ) Q |0}-^ J dt' 



x ( 



+£,j(r < (i / ),t + *>#(*')} 



. 



(521) 



Having substituted (B21) in (B9), we obtain 



(0|W la |0> = ^{0\E a ( Ti ,t)\0) || fjt' 



= exp(iHt'/h){d<p(ri,t + t')/dt') exp(-iHt'/h). 
Then we integrate in parts on t' and obtain 

0*) = "f E<°l^( r *'*)]l°> 



+^(r i (t / ),t + *')^(*')}]|0) 



dt' ( 



jnt'/h 



YyH,ip{v u t + t')\e 



mt'/h 



. 



(518) 



The commutator in the first term is calculated, and for 
transformations of the second term we notice, that 

1 \n-i i + vi 1 ( d^iit) , dtp(r u t) \ 
-[HMrut)} = - 2 \^- B — + - B — Vifi j . (519) 

Having substituted (B19) in (B18), we obtain: 

c^(t) = -^E< i^( r *'*)/ ar -i > + 



2h 



dt' ( 



v-, ,,,Mri(t'),t + t') 

l^Vj*, 2^{Vi0(t )- 



) . (B20) 



2mc ^ \ 



dA a {r t ,t) dA a {r t ,t) 

V i0 g + a,.. ^ 



»/3 



<9r ; 



) . (522) 



This expression coincides with (B8) except for last term 
from the RHS (B22). But it is possible to show, that this 
term is equal 0. Really, 



(1/2) (0 



dA a {Ti,t) dA a (r i7 t) 



dr if3 



= {i/h){0\[H,A a {r i ,t)]\0) = 0, 



(523) 



since the operator 7i has only diagonal matrix elements 
(0|W|0). So, we have checked up, that results (B8) and 
(B9), obtained by different ways, coincide. 

Let us consider a case of free particles, when the con- 
ditions (B7) are carried out. Then C a (t) — 0, the penul- 
timate term from the RHS (B9) also is equal to zero, and 
comparing (B8) and (B4), we find, that for the average 
Lorentz force a ratio owes to be carried out 

^(0|[v< xH(r i ,| a -[H(r i ,t) x Vi ] a |0) /ree = 



tern 
"~2lT 



/" 

J — ( 



dt' ( 
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+Ep(r i (t'),t + t')v if }(t')} 



(524) 



free 



where the index " free" means performance of conditions 
(B7). Let us check up (B24) by direct account, trans- 
forming the RHS. Since Ti.f ree = ( TO /2)X^ w f i v ja com- 
mutes with TLfree, it is possible to write down 



_ JH fm t'/h 



-iHfreet'/H _ 



Vja(t'), 



then a bordering exp(i7i/ I . ee t'/7i) ... — exp(— iHf Tee t' /%) 
from the RHS (B24) leaves. We shall designate the RHS 
(B24) as & a (t), and after calculation of the commutator 
it appears equal 



*«(*) = - 



dr ia 



Vi/3 



Vi/3- 



dr ia 



(525) 



We shall notice, that the equality 



Vip- 



dE a (r u t) , dE a (r h t) 



Or 



ip 



Or 



ip 



-Vi/3 
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Hp 



/ dEp(n,t + f) 

\ dr ia 
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dE a (n,t + t')\ , (dEp{Vi,t + t') 



dE a {r t ,t + t') 



dr- 



ip 



that is equal 

*a(*) = - 



Vip 



free 



(526) 



■ / dt'^2(0\(vi x rot E(ri,t + t')) a 

J — oo 



-(rotE(r i ,t + t / ) x Vi ) a |0). (527) 

Using the Maxwell 

equation rot~E(r,t) = —(l/c)(dH(r,t)/dt) and having 
calculated integral on t', we obtain, that ^ a (t) is equal 
to the LHS (B24), as it was required to prove. 

Let us notice, that in a case E = 0, H = const, 
both sides of (B24) address in 0. It is obvious to 
the RHS, and in the LHS the matrix elements appear 
(0|^i/3|0) = (i/h){0\[H,rip]\Q) = because the operator 
H is diagonal. 
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